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Abstract: A Semi empirical formula for both phase space factor and Nuclear 
Matrix Element (NME) is developed for neutrinoless double beta decay, and the 
formula is used to compute the neutrinoless double beta decay half lives. The 
computed half lives for neutrinoless double beta decay are compared with the 
corresponding experimental values and with those predicted by QRPA model. 
The semi empirical formula predictions are found to be in good agreement with 
experimental data. The semi empirical formula is used to predict neutrinoless 
double beta decay of various isotopes Ca, Ge, Se, Zr, Mo, Pd, Cd, Sn, Te, Xe, Nd 
and Sm that exhibiting single beta decay. As our semi empirical formula predictions 
agree with the experimental data we hope that the present work will be useful for 
the future experiments. 
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1. INTRODUCTION 


Double beta decay is a radioactive decay process where a nucleus releases two 
beta rays as a single process. Here two neutrons in the nucleus are converted in 
to two protons and in the process two electrons and two electron antineutrinos 
are emitted. In order for beta decay to be possible the final nucleus must have 
larger binding energy than the original nucleus. Double beta decay is difficult 
to study in most practically interesting cases, because both beta decay and 
double beta decay are possible, with probability favouring beta decay. The 
double beta decay is usually studied only for beta stable nuclei. Like single 
beta decay, double beta decay does not change the mass number A. More than 
60 naturally occurring isotopes are capable of undergoing double beta decay. 
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Double beta decay is of two types; the two neutrino and neutrinoless 
double beta decay. The two neutrino double beta decay [2(3(2v)] 

(Z, A) —»(Z + 2, A) + e[ + e 2 + v el + v el which involves the 
transformation of two neutrons into two protons conserves not only the electric 
charge but also the lepton number. On the other hand neutrinoless double beta 
decay [2(3(0u)] (Z, A) —> (Z -(- 2, A) + e ] + e 2 violates lepton number 
conservation and is therefore forbidden in the standard electro weak theory. 
According to this theory neutrinos are massless. The observation of neutrino 
mass and oscillation is a clear example of a phenomenon at variance with the 
standard model. 

There are different models for explaining the double beta decay process. 
Among them, two methods are mainly used to calculate Nuclear Matrix 
Elements (NME) for 2(3(Ov) decays. One is the family of Quasi Particle 
Random Phase Approximation (QRPA) [21]. This method has been used by 
different groups and varieties of techniques are employed with results for most 
of the possible emitters [24]. The other method concerned to double beta decay 
process is the Interacting Shell Model (ISM) [8]. It has been shown that as the 
difference in deformation between parent and daughter grows, the NME’s of 
both the neutrinoless and two neutrino mode decreases rapidly. 

The interest in double beta decay spans more than six decades. In 1937 
Racah [16] following the fundamental suggestion of Majorana [11], discussed 
the possibility of a neutrinoless transformation of two neutrons into two 
protons plus two electrons. Even earlier Geoppert-Mayer [28] evaluated the 
decay rate of 2(3(2v) mode and realized that the corresponding half lives could 
exceed 10 20 years. Furry [39] shortly afterwards estimated that 23(0u) should 
be much faster than 2(3(2u) decay. Thus the stage was set for the realization 
that observation of the 2(3(Oy) decay would establish that the neutrino is a 
massive Majorana particle. In 1982 J. Schechter-Valle [23], while considering 
2(3(0u) decay, suggested the existence of Majorana mass of the neutrino in the 
frame work of Gauge theories. In 1984 Fiorini et al [9] introduced a program 
to develop low temperature detectors for 2(3 decay search. Next year Doi et al 
[27] made a fundamental theoretical analysis of 2(3 decay to obtain the main 
formulae for probability of decay, energy and angular electron spectra. In 
1986, using QRPA model Vogel et al [31] got satisfactory agreement between 
theoretical and experimental 2(3(2v) half life values. In the review work of 
Gmez-Cadenas et al [22], the authors provide an answer to the fundamental 
question of whether neutrinos are Dirac or Majorana. 

Neutrinoless double beta decay is of great interest for studying the 
fundamental properties of neutrino beyond the standard electro-weak theory. 
High sensitivity 2(3(0u) studies are the unique and practical ways for studying 
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the Majorana nature of neutrinos, the neutrino mass spectrum, the absolute 
neutrino mass scale, the majorana CP phases and other fundamental properties 
of neutrinos in the foreseeable future. The first experiment to claim 23(0u) is the 
Klapdor, HM [18] experiment done in the year 2001. Numerous experiments 
like COBRA, GERDA etc have been carried out to search neutrinoless double 
beta decay and 48 Ca, 76 Ge, 82 Se, 96 Zr, l00 Mo, 116 Cd, 128 Te, 150 Nd, 238 U are some of 
the isotopes exhibiting neutrinoless double beta decay [4, 5,7,20,32], 

For the double beta decay processes, two crucial ingredients are the phase 
space factors and the Nuclear Matrix Elements (NME). A general theory of 
phase space factors was developed by Doi et al. [25,26] following the previous 
work of Primakoff and Rosen [17], and Konopinski [10]. It was reformulated 
by Tomoda [37] by approximating the electron wave functions at the nuclear 
radius and without inclusion of electron screening. Accurate values of the 
phase space factors are necessary ingredients for theorists, to improve the 
double beta decay lifetime predictions and for experimentalists to plan their 
set-ups [30]. The Nuclear Matrix Element depends on the nuclear structure 
of the nuclei involved in the decay. Frank T Avignone et al [14] have done a 
detailed study on nuclear matrix elements and also pointed out the increasing 
sensitivity of experiments and improvements in nuclear theory make the future 
exciting for this field at the interface of nuclear and particle physics. The 
expression for Nuclear Matrix Element can be written in general as the sum of 
three components [29] as 


M°" = M°J t - 


Sv 


\2 


Sa 


+ M 


0 v 


(1) 


where M^' T , M ( f h ' , Mj V ,aretheGamow-Teller,Fermiandtensorcomponents 
respectively. g A is the axial vector coupling constant and g u is the vector 
coupling constant. 

The present work aims to develop a semi empirical formula for both phase 
space factor and Nuclear Matrix Element for computing the neutrinoless 
double beta decay half life. By using this formula we would like to predict the 
possibility of 2(3(Ou) decay from various isotopes exhibiting single beta decay. 
The details of the semi empirical formula are given in Section 2 and results, 
discussion and conclusion are given in Section 3. 


2. THE SEMI EMPIRICAL FORMULA 

In the standard scenario, when 2(3(Ou) decay process occurs by exchange of 
light Majorana neutrinos between two nucleons inside the nucleus, and in 
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the presence of left handed weak interactions, the life time expression can be 
written as a product of three factors and is given as [3] 

M °" 2 

where G is the phase space factor for this decay mode, (m (/ ) is the effective 
neutrino mass parameter, rn, is the electron mass and M 01 ' are the Nuclear 
Matrix Elements depending on the nuclear structure of the nuclei involved in 
the decay. 

The phase space factor depends on the decay energy Q and nuclear charge 
Z and in the present investigation we studied the dependence of phase space 
factor with ZQ 3 and Z 2 Q 6 for various isotopes undergoing neutrinoless double 
beta decay. From the observed dependence of phase space factor taken from 
ref [34], on ZQ 3 and Z 2 Q 6 , we have developed a semi empirical formula for the 
phase space factor. Using ZQ 3 , Z 2 Q 6 and Z 3 Q 9 as variables, a new formula is 
obtained and is given as, 

= a(ZQ 3 )-b(Z 2 Q 6 ) + c(Z 3 Q 9 ) + d (3) 

The constants are, a = 2.48904 E — 26 , b = 2.2017IE — 38 
c = 9.95199E-51, d= 1.11378E —15 

Due to the two-body nature of the transition operator, the NMEs can also be 
expressed as a sum of product of two-body transition densities (TBTDs) and 
matrix elements of the two-body transition operators for two-particle states. 
We have studied the dependence of nuclear matrix element values taken from 
[38] with Z- 1/3 for various isotopes undergoing neutrinoless double beta 
decay and a new formula is obtained by making least-squares fit to the nuclear 
matrix elements data and is given as, 

M ( 0 °J = aZ~ m + bZ~ m + cZ~ 3/3 + dZ~ 4 ' 3 + eZ ~ 5/3 + / 

The constants are, 

a = —9.49274E + 6 b = 6.65787£ + 7 c = -2.33125E + 8 , 

d = 4.07518£’-t-8 e = -2.84509E + 8 / = 5.40571£ + 5 

The comparison of the computed nuclear matrix elements using the present 
formula with the values of Ref [38] and comparison of computed phase space 
factor with the values of Ref [34] are shown in Table 1. 


7^—1 _/^( 0 ) 

1 1/2 — u iV 
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Table 1. The computed 7j /2 
decay of various isotopes and their comparison' 
and Ref [34] values. 


for neutrino less double beta 
with the experimental, QRPA 


Semi Empirical 
Formula For 
Neutrinoless 
Double Beta Decay 


G ( qJ xio l 5 (y') 


M, 


(0) 

0u 


rp(0v) 
1 1/2 


(yrs) 


Isotope Value 



(KeV) 

Present 

Ref. 

[34] 

Present 

Ref. 

[38] 

Present 

Expt. 

QRPA [38] 

48 Ca 

4272 

24.12 

24.81 

859.27 


5.86E+21 

>5.8E+22[36] 


76 Ge 

2039 

6.44 

2.36 

2.90 

2.48 

1.93E+27 

>1.9E+25[18] 

2.49E+27 

82 Se 

2995 

13.06 

10.16 

3.04 

2.10 

8.63E+26 

>3.6E+23[2] 

8.13E+26 

96 Zr 

3350 

22.59 

20.58 

3.92 

0.40 

3.01E+26 

>9.2E+21[20] 

7.70E+27 

100 Mo 

3034 

16.08 

15.92 

3.48 

1.24 

5.37E+26 

>l.lE+24[20] 

1.45E+27 

nopd 

2018 

7.91 

4.82 

2.57 


2.00E+27 



116 Cd 

2814 

14.73 

16.7 

2.37 

1.31 

1.26E+27 

>1.7E+23[13] 

1.13E+27 

124 Sn 

2287 

10.25 

9.04 

2.38 


1.80E+27 



128 Te 

866 

1.93 

0.59 

2.53 

1.47 

8.47E+27 

>1.5E+24[1] 

2.52E+28 

130 Te 

2527 

12.38 

14.22 

2.53 

1.36 

1.32E+27 

>2.8E+24[7] 

1.15E+27 

136 Xe 

2458 

12.05 

14.58 

2.68 

0.94 

1.21E+27 

>4.5E+23[33] 

2.18E+27 

148 Nd 

1929 

8.54 

10.1 

1.37 


6.54E+27 



150 Nd 

3371 

62.85 

63.03 

1.37 

1.96 

8.88E+26 

>1.8E+22[19] 

1.21E+26 

154 Sm 

1215 

3.62 

3.02 

0.36 


2.24E+29 



160 Gd 

1730 

7.31 

9.56 

3.09 


1.50E+27 

1.30E+21 [12] 


19Sp t 

1047 

3.17 

7.56 

66.26 


7.50E+24 




Table 2. The computed Q values, Phase space factors, nuclear matrix 
elements and the predicted half lives for neutrino less double beta decay of 
various Ca, Ge, Se, Zr, Mo, Pd, Cd, Sn, Te, Xe, Nd, Sm, Gd and Pt isotopes 


Isotope 


Q value 
(KeV) 




T 1/2 (y rs ) 


46 Ca 

988.3 

1.58618E-15 

859.272 

8.92E+22 

47 Ca 

2592.3 

7.53405E-15 

859.272 

1.88E+22 
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Table 2. Continued 


Isotope 

Q value 
(KeV) 

gZ ] (y- 1 ) 

KJ 

T 1/2 (yrs) 

49 Ca 

7269.8 

3.40780E-12 

859.272 

4.15E+19 

50 Ca 

11855.7 

3.44796E-10 

859.272 

4.10E+17 

51 Ca 

13867.8 

1.44900E-09 

859.272 

9.76E+16 

52 Ca 

16955.0 

9.01213E-09 

859.272 

1.57E+16 

75 Ge 

312.6 

1.13809E-15 

2.90139 

1.09E+28 

77 Ge 

3385.6 

1.71285E-14 

2.90139 

7.24E+26 

78 Ge 

5164.1 

5.34962E-13 

2.90139 

2.32E+25 

79 Ge 

6427.6 

4.72938E-12 

2.90139 

2.62E+24 

80 Ge 

8244.9 

5.07807E-11 

2.90139 

2.44E+23 

81 Ge 

10089.5 

3.30371E-10 

2.90139 

3.76E+22 

83 Se 

4640.7 

2.21947E-13 

3.04409 

5.08E+25 

84 Se 

6479.0 

6.21762E-12 

3.04409 

1.81E+24 

85 Se 

9052.0 

1.46232E-10 

3.04409 

7.71E+22 

86 Se 

12724.6 

3.31456E-09 

3.04409 

3.40E+21 

87 Se 

14129.4 

8.57934E-09 

3.04409 

1.31E+21 

88 Se 

15812.0 

2.37729E-08 

3.04409 

4.74E+20 

94 Zr 

1142.9 

2.52372E-15 

3.91925 

2.69E+27 

95 Zr 

2049.7 

7.48180E-15 

3.91925 

9.09E+26 

97 Zr 

4593.8 

3.46849E-13 

3.91925 

1.96E+25 

98 Zr 

6824.7 

1.72152E-11 

3.91925 

3.95E+23 

"Zr 

8197.8 

9.63615E-11 

3.91925 

7.06E+22 

100 Zr 

9584.0 

4.08099E-10 

3.91925 

1.67E+22 

101 Zr 

10051.0 

6.31459E-10 

3.91925 

1.08E+22 
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Table 2. Continued 


Isotope 

Q value 
(KeV) 

G" (y-'> 

KJ 

T i/2 (y rs ) 

10 “Zr 

11817.0 

2.76762E-09 

3.91925 

2.46E+21 

103 Zr 

12480.0 

4.54642E-09 

3.91925 

1.50E+21 

98 Mo 

112.3 

1.11526E-15 

3.47678 

7.75E+27 

"Mo 

1651.2 

5.10021E-15 

3.47678 

1.69E+27 

101 Mo 

4438.7 

2.88638E-13 

3.47678 

2.99E+25 

102 Mo 

5541.0 

2.68528E-12 

3.47678 

3.22E+24 

103 Mo 

6408.8 

1.10330E-11 

3.47678 

7.83E+23 

104 Mo 

7759.0 

6.71620E-11 

3.47678 

1.29E+23 

105 Mo 

8588.0 

1.72446E-10 

3.47678 

5.01E+22 

106 Mo 

10067.0 

7.43636E-10 

3.47678 

1.16E+22 

107 Mo 

11430.0 

2.37009E-09 

3.47678 

3.65E+21 

109p d 

901.0 

1.92669E-15 

2.56879 

8.22E+27 

11 ‘Pd 

3253.5 

2.48567E-14 

2.56879 

6.37E+26 

“2p d 

4244.5 

2.49383E-13 

2.56879 

6.35E+25 

113 Pd 

5359.3 

2.60657E-12 

2.56879 

6.07E+24 

n-tpd 

6523.9 

1.74643E-11 

2.56879 

9.06E+23 

115 Pd 

7690.5 

8.20373E-11 

2.56879 

1.93E+23 

“ 6 Pd 

8759.0 

2.73683E-10 

2.56879 

5.78E+22 

117 Pd 

9895.0 

8.38279E-10 

2.56879 

1.89E+22 

118 Pd 

11239.0 

2.67934E-09 

2.56879 

5.91E+21 

ll4 Cd 

540.1 

1.30076E-15 

2.37189 

1.43E+28 

115 Cd 

1945.5 

7.60027E-15 

2.37189 

2.44E+27 

“ 7 Cd 

3975.0 

1.46594E-14 

2.37189 

1.25E+26 

118 Cd 

4947.1 

1.48733E-13 

2.37189 

1.37E+25 

ll9 Cd 

6158.4 

1.35555E-12 

2.37189 

1.61E+24 


Semi Empirical 
Formula For 
Neutrinoless 
Double Beta Decay 


171 










Raj an, MKP 
Biju, RK 
Santhosh, KP 


Table 2. Continued 


Isotope 

Q value 
(KeV) 

(y- 1 ) 

KJ 

T 1/2 (yrs) 

120 Cd 

7131.1 

1.15360E-11 

2.37189 

4.01E+23 

121 Cd 

8144.1 

4.62494E-11 

2.37189 

1.17E+23 

122 Cd 

9215.9 

1.59314E-10 

2.37189 

3.73E+22 

123 Cd 

10510.5 

4.97662E-10 

2.37189 

1.12E+22 

124 Cd 

11526.8 

1.65584E-09 

2.37189 

4.84E+21 

125 Cd 

12538.5 

3.83662E-09 

2.37189 

2.25E+21 

126 Cd 

13690.0 

8.23556E-09 

2.37189 

1.02E+21 

127 Cd 

14979.0 

1.82585E-08 

2.37189 

4.50E+20 

128 Cd 

16045.0 

4.12121E-08 

2.37189 

2.42E+20 

129 Cd 

17394.0 

7.67089E-08 

2.37189 

1.17E+20 

130 Cd 

20722.0 

1.59010E-07 

2.37189 

2.41E+19 

131 Cd 

22044.0 

7.71382E-07 

2.37189 

1.38E+19 

132 Cd 

25834.0 

1.34703E-06 

2.37189 

3.30E+18 

122 Sn 

368.1 

1.17572E-15 

2.38151 

1.57E+28 

123 Sn 

1351.4 

3.86873E-15 

2.38151 

4.76E+27 

125 Sn 

3123.7 

2.31358E-14 

2.38151 

7.96E+26 

126 Sn 

4044.6 

2.02822E-13 

2.38151 

9.08E+25 

127 Sn 

4782.1 

1.10593E-12 

2.38151 

1.67E+25 

128 Sn 

5657.1 

5.80415E-12 

2.38151 

3.17E+24 

129 Sn 

6409.2 

1.92150E-11 

2.38151 

9.58E+23 

130 Sn 

5059.4 

1.94125E-12 

2.38151 

9.49E+24 

131 Sn 

7895.5 

1.35612E-10 

2.38151 

1.36E+23 

132 Sn 

8628.0 

3.07712E-10 

2.38151 

5.98E+22 

133 Sn 

11995.0 

6.23293E-09 

2.38151 

2.95E+21 

134 Sn 

15759.0 

7.37341E-08 

2.38151 

2.50E+20 

135 Sn 

17030.0 

1.48546E-07 

2.38151 

1.24E+20 

136 Sn 

17930.0 

2.36434E-07 

2.38151 

7.79E+19 

137 Sn 

19250.0 

4.48748E-07 

2.38151 

4.10E+19 
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Table 2. Continued 


Isotope 

Q value 
(KeV) 

C <r‘) 

<? 

T i/2 (y rs ) 

127 Te 

39.9 

1.11386E-15 

2.52753 

1.47E+28 

129 Te 

1694.2 

6.16089E-15 

2.52753 

2.65E+27 

131 Te 

3205.7 

2.91713E-14 

2.52753 

5.60E+26 

132 Te 

4098.5 

2.64656E-13 

2.52753 

6.18E+25 

133 Te 

4698.6 

1.05662E-12 

2.52753 

1.55E+25 

134 Te 

5565.5 

5.62396E-12 

2.52753 

2.91E+24 

135 Te 

8587.0 

3.32169E-10 

2.52753 

4.92E+22 

136 Te 

11995.0 

7.0181 IE-09 

2.52753 

2.33E+21 

137 Te 

12819.0 

1.28185E-08 

2.52753 

1.28E+21 

138 Te 

14220.0 

3.27787E-08 

2.52753 

4.99E+20 

139 Te 

14844.0 

4.83315E-08 

2.52753 

3.38E+20 

140 Te 

16030.0 

9.67914E-08 

2.52753 

1.69E+20 

141 Te 

16770.0 

1.45479E-07 

2.52753 

1.12E+20 

142 Te 

18050.0 

2.82536E-07 

2.52753 

5.79E+19 

134 Xe 

825.4 

1.84958E-15 

2.68156 

7.85E+27 

l35 Xe 

1433.5 

4.55604E-15 

2.68156 

3.19E+27 

137 Xe 

5342.2 

4.26727E-12 

2.68156 

3.40E+24 

139 Xe 

9269.7 

7.52261E-10 

2.68156 

1.93E+22 

140 Xe 

10281.0 

1.93663E-09 

2.68156 

7.50E+21 

141 Xe 

11396.0 

4.94136E-09 

2.68156 

2.94E+21 

142 Xe 

12343.0 

1.01957E-08 

2.68156 

1.42E+21 

143 Xe 

13486.0 

2.27398E-08 

2.68156 

6.39E+20 

144 Xe 

14489.0 

4.35099E-08 

2.68156 

3.34E+20 

145 Xe 

15310.0 

7.15996E-08 

2.68156 

2.03E+20 

146 Xe 

16330.0 

1.28198E-07 

2.68156 

1.13E+20 

147 Xe 

17340.0 

2.20354E-07 

2.68156 

6.59E+19 

146 Nd 

70.9 

1.11431E-15 

1.36803 

5.01E+28 
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Table 2. Continued 


Isotope 

Q value 
(KeV) 

(y- 1 ) 

Ki 

T 1/2 (yrs) 

147 Nd 

1120.2 

3.0624IE-15 

1.36803 

1.82E+28 

149 Nd 

2761.0 

1.74770E-14 

1.36803 

3.19E+27 

151 Nd 

3629.5 

1.26279E-13 

1.36803 

4.42E+26 

l52 Nd 

4610.8 

1.41057E-12 

1.36803 

3.96E+25 

l53 Nd 

5216.8 

4.76726E-12 

1.36803 

1.17E+25 

154 Nd 

6771.6 

5.71782E-11 

1.36803 

9.76E+23 

155 Nd 

7727.2 

1.94953E-10 

1.36803 

2.86E+23 

156 Nd 

8840.0 

6.71862E-10 

1.36803 

8.31E+22 

157 Nd 

9940.0 

1.96131E-09 

1.36803 

2.85E+22 

158 Nd 

10810.0 

4.20847E-09 

1.36803 

1.33E+22 

159 Nd 

11990.0 

1.07758E-08 

1.36803 

5.18E+21 

160 Nd 

13000.0 

2.24165E-08 

1.36803 

2.49E+21 

161 Nd 

14020.0 

4.43901E-08 

1.36803 

1.26E+21 

l53 Sm 

324.0 

1.16617E-15 

0.35918 

6.94E+29 

156 Sm 

3172.2 

4.12843E-14 

0.35918 

1.96E+28 

157 Sm 

4100.7 

4.82787E-13 

0.35918 

1.68E+27 

l58 Sm 

5486.8 

8.63633E-12 

0.35918 

9.37E+25 

159 Sm 

6358.5 

3.51019E-11 

0.35918 

2.31E+25 

160 Sm 

7528.6 

1.69543E-10 

0.35918 

4.78E+24 

161 Sm 

8532.7 

5.36970E-10 

0.35918 

1.51E+24 

162 Sm 

9537.0 

1.48573E-09 

0.35918 

5.45E+23 

163 Sm 

10590.0 

3.85573E-09 

0.35918 

2.10E+23 

164 Sm 

11570.0 

8.61159E-09 

0.35918 

9.40E+22 

165 Sm 

12670.0 

1.96091E-08 

0.35918 

4.13E+22 

159 Gd 

605.0 

1.46214E-15 

3.090977 

7.48E+27 

161 Gd 

2548.4 

1.46029E-14 

3.090977 

7.49E+26 

162 Gd 

3899.8 

3.22654E-13 

3.090977 

3.39E+25 

163 Gd 

4896.5 

3.16629E-12 

3.090977 

3.45E+24 

164 Gd 

6223.3 

3.16750E-11 

3.090977 

3.45E+23 
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1S5 Gd 

7147.9 

1.15629E-10 

3.090977 

9.45E+22 

Semi Empirical 
Formula For 






Neutrinoless 

166 Gd 

8190.1 

4.06225E-10 

3.090977 

2.69E+22 

Double Beta Decay 

167 Gd 

9240.0 

1.22598E-09 

3.090977 

8.92E+21 


168 Gd 

10460.0 

3.79423E-09 

3.090977 

2.88E+21 


169 Gd 

11700.0 

1.04894E-08 

3.090977 

1.04E+21 


197 Pt 

118.6 

1.11702E-15 

66.26523 

2.13E+25 



199pt 

2155.1 

1.18621E-14 

66.26523 

2.01E+24 

200pt 

2901.1 

3.74061E-14 

66.26523 

6.36E+23 

201 Pt 

3923.1 

6.69574E-13 

66.26523 

3.55E+22 

202p t 

4745.9 

4.44652E-12 

66.26523 

5.35E+21 


3. RESULTS, DISCUSSION AND CONCLUSION 

The Q value for double beta decay of mother nuclide with mass m m to the 
daughter nuclide with mass m is given by the mass difference Q = m m — m d 
[6] which in turn can be written as a function of the ratio of cyclotron frequencies 
of the ions in a Penning trap [35], r = v d / v and the electron mass m e ; 


Q = tn m - m, = 


1 

V 

m 


(m d -m e ) 


(5) 


In the present work Q values are computed using the experimental binding 
energies of Audi and Wapstra [15]. The present empirical formula is applied 
for all the observed neutrinoless double beta decay isotopes. Column 7 of 
Table 1 represents the computed half-lives for neutrinoless double beta decay 
of various isotopes and is compared with the experimental values given in 
column 8 and QRPA values [38]in column 9. It is found from the table that 
our formula predictions are in good agreement with the experimental values 
and the QRPA values. The value of <m> is taken as 50MeV and is obtained 
from Rodin et al [38]. 

We have applied the present formula for computing the phase space factor, 
Nuclear Matrix Element and half lives for various isotopes that exhibiting 
single beta decay. Tables 2 represents the computed Q values, Phase space 
factors, Nuclear Matrix Elements and half lives for neutrinoless double beta 
decay of various Ca, Ge, Se, Zr, Mo, Pd, Cd, Sn, Te, Xe, Nd and Sm isotopes. 
As our semi empirical formula prediction agree with the experimental data we 
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hope that our prediction on neutrinoless double beta decay of various Ca, Ge, 

Se, Zr, Mo, Pd, Cd, Sn, Te, Xe, Nd and Sm isotopes will be a guide for future 

experiments. 
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